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The anomalous response of glasses in the echo amplitude experiment is explained in the presence 
of a magnetic field. We have considered the low energy excitations in terms of an effective two level 
system. The effective model is constructed on the flip-flop configuration of two interacting two level 
systems. The magnetic field affects the tunneling amplitude through the Aharanov-Bohm effect. The 
effective model has a lower scale of energy in addition to the new distribution of tunneling parameters 
which depend on the interaction. We are able to explain some features of echo amplitude versus a 
magnetic field, namely, the dephasing effect at low magnetic fields, dependence on the strength of 
the electric field, pulse separation effect and the influence of temperature. However this model fails 
to explain the isotope effects which essentially can be explained by the nuclear quadrupole moment. 
We will finally discuss the features of our results. 
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I. INTRODUCTION 

Unusual behavior of the heat capacity and thermal 
conductivity of glasses at very low temperatures^ can 
be governed by a broad distribution of low-energy ex- 
citations. These could be described as two level tun- 
neling systems whose energies and tunneling amplitudes 
are randomly distributed^. The origin of these excita- 
tions is assumed to come from atoms or group of atoms 
that can be arranged in two energetically close configura- 
tions. Such a configuration can be a cluster of few hun- 
dred atoms or molecules^. Although the nature of the 
tunneling entities is obscure, many theories are based on 
such a tunneling model to describe the low temperature 
properties of glasses. A two level system (TLS) is char- 
acterized by two parameters of a double well potential, 
the asymmetry energy A, and the tunneling amplitude 
Ao . The Hamiltonian of a single TLS in the pseudo-spin 
operators is given by 

H = -(Aa z + A a x ), (1) 

As a consequence of the irregular structure of glasses, 
these two parameters are widely distributed and are in- 
dependent of each other with a uniform distribution of 

/«(A,A ) = A (2 ) 

where fo is a normalization factor. While many exper- 
iments can be explained satisfactorily by this model it 
fails to explain some results at very low temperatures, 
particularly in the presence of a magnetic field^. 

It was a general belief that the dielectric properties 
of nonmagnetic insulating glasses are independent of the 
magnetic field. Some experiments during the last few 
years opened the question of magnetic field dependence. 
The dielectric constant of several multicomponent glasses 



like BaO — AI2O3 — S1O2 shows a surprising dependence 
on a weak magnetic field (~ 20/zT) at very low tempera- 
tures (~ 2mK)&L. Increasing the magnetic field the de- 
pendence shows an oscillatory behavior with a period of a 
few Teslas. Moreover the amount of magnetic impurities 
present in the nonmagnetic glasses seems to be irrelevant 
and the observed magnetic field effect is strongly influ- 
enced by the applied electric fields. On the other hand, 
recent polarization echo experiments show a strong and 
nonmonotonic dependence of echo amplitude on the ap- 
plied magnetic field^S. These experiments give us an 
insight into the coupling of the tunneling centers to the 
magnetic field. Namely, the echo amplitude is almost 
independent of the frequency of the electric field which 
means that tunneling systems with very different energy 
splittings behave the same waj*ii. Moreover the effect is 
strongly dependent on the delay time between two pulses 
and also on the magnitude of applied electric field. A sur- 
prising outcome of these experiments is a novel isotope 
effect observed in different glasses^. The recent results 
show the important influence of nuclear quadrupole mo- 
ments on the observed magnetic field dependence. 

The first theoretical attempt to explain the mag- 
netic field dependence was based on the Aharanov-Bohm 
phased. In this approach the dielectric permittivity is 
given by a generalized TLS where the magnetic field 
enters in the quantum mechanical phase of a tunneling 
charged particle in a Mexican hat potential. It was as- 
sumed that such a generalized TLS is the result of collec- 
tive effects between many tunneling centers in a meso- 
scopic size which carry a big amount of charge. The 
dipole-dipole interaction between tunneling centers is re- 
sponsible for the mentioned collective (coherent) mo- 
tion in a weakly interacting regimei^. An evidence for 
the Aharanov-Bohm effect can be found in the metal- 
lic glasses^. Another approach which is also based on 
the Aharanov-Bohm phase explains the phenomena of 
dielectric response in terms of coupled pairs of tunneling 



2 



centers ™. The coupling of two TLSs makes the tunneling 
path curved and possessing a flux dependent phase. The 
most recent theory is based on the nuclear quadrupole 
moment of tunneling entitiesiilS. In this model the in- 
teraction of nuclear quadrupole moment with the gradi- 
ent of the electric field defines the new scale of energy at 
small magnetic fields. The different orientations of the 
nuclear quadrupole moment in each well lead to differ- 
ent level splittings which affect the dephasing mechanism 
and finally the echo amplitude. This model explains the 
echo experiments well, however, it fails to get the right 
order of magnitude for the magnetic field effect of di- 
electric response^. It is claimed that the observed effect 
in dielectric response is not a static phenomenon but of 
dynamical origin^. 

The shortcoming of the isolated TLS model to explain 
physical phenomena at low temperatures is usually re- 
lated to the interaction between tunneling centersSi 2 ^. 
Especially, interaction plays an essential role in deter- 
mining the relaxation rates at low temperatures 2 - 3 which 
influences the echo experiment response. In this most re- 
cent work it has been shown how the dipole-dipole inter- 
action can lead to a relaxation rate faster than phononic 
counterpart at very low temperatures. The main idea 
is based on the derealization process of the effective 
TLSs which are the result of two interacting tunneling 
centers^S^. This is in agreement with the experimental 
results on the magnetic field dependence of the dielectric 
constant and the effective model proposed in Ref . [2^ . 
Here, the magnetic field dependence of the dielectric con- 
stant is related to the field dependence of the tunneling 
amplitude cutoff which is based on the theory developed 
earlier 2 ' . In this theory the localization length increases 
by adding the magnetic field which favors the derealiza- 
tion phenomenon. 

Accepting the fact that interaction plays an essential 
role in the low temperature physics of glasses we are go- 
ing to explain the recent polarization echo experiments in 
the presence of a magnetic field. Our appro ach is based 
on the main idea presented in Refs. 13 23 26]. The inter- 
action between two tunneling centers is defined in terms 
of flip-flop states which are considered as a new (effective) 
TLS. The magnetic field adds the Aharanov-Bohm phase 
to the tunneling process. The effective TLS responds to 
the electric field of two pulses by means of renormalized 
parameters and distribution. We are able to get most of 
the features of echo experiments, namely, strong depen- 
dence of the echo amplitude on low magnetic fields, pulse 
separation, temperature and electric field effects. Our re- 
sults show how one can explain the echo experiments by 
using the idea of Aharanov-Bohm phase. 

In the next section we define our effective model and 
its renormalized parameters. The response of an effective 
TLS to a two pulse echo is studied in Sec lIIII Then, in 
Sec II VI we will explain the different aspects of our results 
to show their correspondence with experimental data. Fi- 
nally we summerize our results and discuss the features 
of our approach. 
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FIG. 1: Flip-flop configuration as an effective TLS. 

II. EFFECTIVE TLS 

Suppose that there are two interacting tunneling sys- 
tems. The interaction between the TLSs can be consid- 
ered by the Ising Hamiltonian 
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U(Rij) 



1,3 



(3) 



where Rij is the distance between two TLSs and Uq is 
the characteristic coupling constant. 

Burin et ali^i proposed that TLS-TLS coupling in 
glasses is dominated at low temperatures by a flip-flop 
configuration. While the characteristic energy of a sin- 
gle TLS (E or E', see Fig.Q) has the same order of 
magnitude as temperature (T), the resulting spectrum 
of two TLSs has two energy scales. The flip-flop states 
define a scale of energy which is much less than T and 
the other remaining two represent a scale as before, i.e. 
T. The flip-flop (pair) states correspond to the situation 
where one TLS is in the ground state and the other in the 
excited state. On the other hand the remaining part of 
Hilbert space shows both excited or both grounded TLSs. 
Concerning the relaxation of the systems which leads fi- 
nally to the dephasing mechanism we only consider the 
flip-flop states. Taking into account the interaction be- 
tween tunneling centers the Hamiltonian in the flip-flop 
subspace is 



R p = -(A p cr z + A 0p cr x ) 



where 



u 



AqiAc 
EE' 



A P =\E- E' 



(4) 



The distribution function for the parameters Ao p and 
A p is defined as 

/( 2 )(A p ,A 0p )= <5(A P -\E-E'\)x 

^ - i?3— EtfH > " (5) 
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Therefore by averaging over the distribution parameters 
given by Eq. (J2J of the original TLS and ensemble aver- 
aging (A p < E ~ T) we get 



/( 2 )(A p ,A 0p )cx-g-. 



(6) 



Now instead of a pair of original TLSs we have an ef- 
fective TLS (ETLS) which behaves like a single TLS with 
the new distribution function, Eq.©. Such an effective 
TLS is aimed to explain the response of a glass to an 
electric field of the echo experiment at low temperatures 
when a magnetic field is present. The magnetic field de- 
pendence enters via the dependence of the original tun- 
neling amplitude of pairs (Aoi(</>) = Aoi cos(7r^-); i — 
1,2) as proposed in Ref . [l3| . By a simple calculation 
one can show that the magnetic field dependence of 
Ao p is similar to the original TLS one's, i.e. Ao p (<fr) « 
Ao p cos(7r^-), where <fio is the renormalized quantum of 
flux (h/q). The Planck constant is denoted by h and the 
effective charge is q. We should mention that this model 
is valid while the temperature is not zero. At the zero 
temperature the distribution function of ETLSs, Eq.©, 
is zero which means no accessible flip-flop states. In other 
words, at T — the system does not respond to the elec- 
tric field of the echo experiments and will remain in the 
ground state. 



III. THE RESPONSE OF AN EFFECTIVE TLS 
TO A TWO PULSE ECHO 

The wave function of an ETLS in an external time 
dependent electric field is a linear combination of the 
wave function \ipi) and \tp 2 ) for each level. 



| ip) = ai(t)e~ 



ipi) + a 2 (t)e 



-iE 2 t 



(7) 
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FIG. 2: The echo amplitude versus the magnetic flux ratio 
for relaxation time to = 2fis, by averaging over the effective 
distribution function for < % < lQ s Hz and 10 9 < %^ < 



10 9 + 10 6 JLz. 
15mK). 
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the ETLS, respectively. We can easily find that at t = t\ 
the coefficients are given by, 



oi (n) 



a 2 {Ti) 
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7i + 72 
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(9) 



|o x (t) \ 2 + \a 2 (t) | 2 =1. 



At any time (t) during the first pulse of the electric field 
(£(t) — ^ sin cut), the time variation of the probability 
amplitudes a\(t) and a 2 (t) is governed by the Schrodinger 
equation. In the resonant condition (| hui — E |<C fuv) it 
leads to the following equations^, imposing A p -C Aop^, 

dai x iSt m 
at 



da 2 



=± = -* e- wt ai(*), 



(8) 



Po-g^ in which p 



where 5 = u — ¥■ , So = and pi 2 
is the dipole moment of the ETLS. E = E 2 — E\ where E\ 
and E 2 are the energy of ground and the exited states of 



where 7l = S W^K and 72 = z£±Vp±M. The 
initial condition is imposed by the Boltzmann weight, 
i.e. |ai(0)| 2 = ^e"^ 1 and |a 2 (0)| 2 = ^e~ 0E \ where 



-/3-Ei 



-0E 2 



is the partition function, (3 



K B T 



and Kb is the Boltzmann constant. 

The system relaxes for a time duration of to before the 
second pulse arrives. The second pulse of the electric field 
tries to recover the initial configuration during t 2 = 2t\ 
seconds. Finally the measurement on the amplitude of 
the electric dipole is done at t' = to after switching off 
the second pulse. This can be expressed by the following 
equation, 



A(2t ) = < i/>(2to) I p I V(2*o) > 
~ 2Re[al(T 2 )a 2 (T 2 )pi 2 e~' 



(10) 



where p = pocr z and the details of the other coefficients 
can be found in the appendix. 
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FIG. 3: The echo amplitude versus the magnetic flux ratio 
for different electric field amplitudes, by averaging over the 
effective distribution function for < ^ < 10 s Hz and 10 9 < 
|f < 10 9 + 10 6 Hz, (lj = 2 x 10 9 Hz, to = 2fis and T = 
15mK). The vertical axis has been shifted. 



FIG. 4: The echo amplitude versus electric field amplitude, 
for different values of the magnetic field flux, by averaging 
over the effective distribution function for < A < W 7 Hz 
and 10 9 < |jf < 10 9 + 10 6 Hz, (d = 2 x 10 9 //z, to = 1/xs and 
T = 15mK). 



As mentioned before we will consider the effect of the 
magnetic field via a Mexican hat model. In this model 
for describing the non-monotonic magnetic field depen- 
dence of the dielectric susceptibility of multi-component 
glasses, Kettemann et al^ investigated the properties 
of tunneling systems exhibiting the peculiarity that the 
tunneling particle can move along different paths to go 
from one potential minimum to the other, thus forming 
a closed tunneling loop. As a result in the three dimen- 
sional Mexican hat model, using the Aharanov-Bohm ef- 
fect and the assumption that the tunneling through both 
paths occurs with equal probability, the tunnel splitting 
becomes a periodic function of the magnetic flux thread- 
ing through the loop, 



A 0p O) = A pcos(7r— ), 



(11) 



where </3 is the total magnetic flux passing through a 
closed tunneling path, cf>o is the renormalized quantum 



flux defined by 



■| in the previous section. 



when the amplitude of the driving electric field is given 

bysa 



2h ~ 4 



(12) 



and at T = lQmK the tunneling matrix element is ^ ~ 
10 9 Hz. If we consider the typical value for £ Q = 1 — 
and ti = 100ns then £^ - 10 6 - W 7 Hz. To compare 
our results with recent experiments2*i2*ii*i^ the following 
frequency ^ ~ 10 8 — 10 9 Hz and pulse separation t = 
t a ~ 10~ 6 s are assumed. 

We should mention that the result of Ea. I|10fl is valid 
for low magnetic field amplitudes. At high magnetic 
fields the effective tunneling amplitude A 0p cos(7T^) be- 
comes small and the resonance condition (|<5| <^ ui) is 
not fulfilled. Therefore, the resonance condition (|w — 

^/Al+A 0p (4>) 2 | < ^ ig satisfied wnenever 1 _ C0S(7T^) < 

1, having assumed that ~ u>. In other words the 
resonance condition is guaranteed for (|7r-?-| -C I). 



IV. RESULTS ON ECHO AMPLITUDE 

In this section we are going to present our results on the 
amplitude of the echo response (Ea. HlOjl ) as a function of 
the applied magnetic and electric field in addition to the 
pulse separation time (to) and temperature. 

In the absence of a magnetic field the echo amplitude 
measured for small pulse separations in an ideal double- 
pulse experiment, ignoring relaxation, has a maximum 



A. Magnetic field dependence 

We have plotted in Fig.@, the square of the normal- 
ized echo amplitude (\A(B)/A(B = 0)| 2 ) versus the mag- 
netic field. The magnetic field (B) in the horizontal axis 
is presented as flux ((f) = SB, S is the effective area of 
ETLS) divided by the quantum flux ((/>o) to simplify data 
analysis. To arrive at this plot we have used Ea. Hl(J|) 
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FIG. 5: The echo amplitude versus the magnetic flux ratio for 
different delay times t — 1/is and t = 1.5/is, by averaging 
over the effective distribution function for < £ < IQpHz 
and 10 9 < |f < 10 9 + 10 6 #2, (w = 2 x 10 9 #z, ^ = 
7 x 10 6 #z and T = 15mK). 



FIG. 6: The echo amplitude versus the magnetic flux ratio 
for different delay times, by averaging over the effective dis- 
tribution function for < f < W a Hz and 10 9 < §f < 
10 9 + 10 6 Hz, ((j = 2 x W 9 Hz, = 7 x 10 6 #2 and 

T = I5mi\). The vertical axis has been shifted. 



with to = 2/xs. We have also averaged over the tun- 
neling parameters (A p , A 0p ). Note that, in this process 
we have assumed A p < Ao p to have reasonable flip-flop 
configurationsSiSi. Moreover the resonance condition of 
the ETLS limits the range of Ao p to be very close to the 
frequency of the electric field 10 9 < ^ < 10 9 + 10 6 Hz. 
The distribution of the ETLS parameters is given by 
Eq.@. The echo amplitude has a maximum at B = 
as presented in Fig.|j2J). It then decreases to a minimum 
before rising up. This is in agreement w ith th e experi- 
mental observation presented in Ref.|fl ll(Mllll2j . 



B. Dependence on the electric field 

We have examined the dependence of our results on 
the amplitude of the electric field. The echo amplitude 
versus the magnetic field for different strengths of the 
electric field is plotted in Fig.©. The data for each value 
of the electric field have been shifted by a constant for 
clarity. We will see that the height of the central peak 
is enhanced by increasing the electric field. This feature 
also coincides with the experimental data on the defected 
crystal KBr : CN (Fig.l in Ref.0). 

We have then plotted the echo amplitude versus the 
amplitude of the electric field in Fig. Q for different mag- 
netic fields. The echo amplitude rises and has a maxi- 
mum for an electric field around 700V/m, it then de- 
creases for larger fields. The general feature is in agree- 
ment with the experiments on multicomponent glasses 
(Fig.2 in Ref.@). 



C. Effect of pulse separation 

We have presented in Fig.© the same quantity as 
Fig-© but for different pulse separations (to = 1, 1.5/xs). 
Although the global behavior is similar to Fig.© there 
is a slight difference at small magnetic fields. We have 
observed that for short pulse separations, to = 1/is, the 
echo amplitude gets a minimum at B = 0. The effect is 
similar to the experiments on crystals with point defects 
(Fig. 2 of Ref.[ji|). For short pulses the magnetic field 
does not have a dephasing effect. However the central 
peak broadening does not accompany the decreasing of 
pulse separation as in Fig. 3 of Ref . |10| . In contrast, the 
broadening of the central peak decreases with to- This 
can be seen in Fig. 10 where we have scaled the horizon- 
tal axis by a factor of to. The width of the central peak 
diminishes by applying shorter pulse separations. 



D. Temperature dependence 

The temperature dependence of our results enters via 
the Boltzmann weight as the initial probability of the 
states in ETLS. This dependence is plotted in Fig.JTJ). 
We have plotted the amplitude of the echo response ver- 
sus magnetic field for different temperatures. To have 
a clear plot for different curves we have added a con- 
stant value to each curve. Thus they are well represented 
instead of falling onto each other. We have seen that 
the intensity of the central peak is reduced by increasing 
the temperature and becomes almost flat at high tern- 
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FIG. 7: The echo amplitude versus the magnetic flux ratio 
for Different temperatures, by averaging over effective dis- 
tribution function for < £ < IVPHz and 10 9 < §f < 
10 9 + 10 6 //z, (w = 2xl0 9 Hz, t = 2/isand = TxltfHz). 
The vertical axis has been shifted. 



peratures. Apart from that the general behavior is un- 
changed. This can also be compared to the experimental 
data on the multicomponent glasses (Fig. 8 of Ref . . 
Although the profile of the echo amplitude versus mag- 
netic field does not fit completely, the reduction of central 
peak is in agreement. 



V. SUMMARY AND DISCUSSION 

We have presented a model to explain the anomalous 
behavior of the echo response in glasses at low temper- 
atures. This model is essentially based on an effective 
TLS which is the result of interaction between the orig- 
inal TLSs in glasses or defected crystals. Most of the 
experimental features at low magnetic fields can be ex- 
plained by our approach. In this respect, the dephas- 
ing effect which reduces the echo amplitude at low fields 
(Fig.®), the dependence on the strength of the electric 
field (Figs. (PI 0J)) and also the effect of various temper- 
atures (Fig. 10) have been addressed. However, in the 
response of different pulse separations (to) although it 
shows the broadening of the central peak by reduction of 
to the product of to4> does not seem to be universal which 
deviates from the experimental data. 

As we have already mentioned our model is valid for 
the low magnetic field regime. At high magnetic fields 
the resonance condition is lost, thus the flip-flop basis of 
ETLS does not contribute effectively. Because, if a TLS 
is initially in the ground state, the probability to find the 



A 2 + A (# 



•sin 2 ( 



VA 2 + A o (0) 2 -^ 

2h 



(13) 



For high magnetic fields this probability vanishes since 
Ao(0) become very small. So the particles will be frozen 
in the lowest well and do not respond to the external 
electric field. Thus we expect the echo amplitude for 
high magnetic fields to be constant. 

The other aspect of echo response which can be ex- 
plained by this model is the response to the frequency 
of the electric field. It has been shown that the echo 
amplitude of BK7 is almost independent of electric field 
frequency (see Fig. (7) in Ref.jTJ]). A raw description 
is given by taking into account the interaction between 
ETLS defined in Secjnj Thus the resonance condition 
can take place for different frequencies and the final re- 
sponse will be the same. Note that the distribution of 
the effective ETLS is given by Eq.© which is given in 
the appendix. 

However, there is a debate on the interpretation of the 
results which comes from a Mexican hat model. Let us 
reduce the discussion to two points. The first point is 
related to the route where the magnetic field enters the 
problem. In our model (Mexican hat type) the magnetic 
field effect appears in the tunneling amplitude (Eq.jnj) 
via the Aharanov-Bohm phenomenon. This by itself can 
not be ruled out. Recently, Cochec, et &l~ reported ex- 
periments on the dielectric constant of a structural glass 
and explained their results by considering the variation 
of the tunneling amplitude versus magnetic field. This is 
based on a microscopic theory of transport in disordered 
networks which can be interpreted as the tunneling be- 
tween two wells. The effect of the magnetic field due to 
quantum phases influences the tunneling amplitude^. 

The second point is related to the interpretation of a 
large effective charge q, (or small quantum flux fa) which 
scales the magnetic field to the right value of the experi- 
mental data. This is the most controversial point and we 
are going to add some comments in this respect. Firstly, 
the large value of q has been taken as an indication that 
there is a coherent motion of tunneling particles taking 
place as the result of their mutual interactioniii4. As 
far as the low temperature physics of glasses is concerned 
the interaction between TLSs is an important issue which 
has been recently addressed as a dephasing mechanism 
in the relaxation of cold glasses^. The relaxation can be 
enhanced by interaction through flip-flop configurations 
leading to delocalized state. Thus the interaction should 
also have a considerable effect in the dephasing mecha- 
nism in the presence of a magnetic field. And this is taken 
into account in the Mexican hat type model effectively. 
Secondly, the nature of a tunneling system (TLS) is not 
clearly known. The general belief is that it represents the 
low energy excitations of glasses at low temperatures. 
It has been recently argued that such types of excita- 
tion come from two energetically close configurations of 
a cluster of atoms or molecules^S. The cluster contains 
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roughly N — 200 atoms or molecules contributing in a 
collective tunneling. Thus the flux passing through such 
a cluster is much larger than the case in a single atom. 
This by itself can improve the necessary flux (<fi) in the 
Mexican hat model by two orders of magnitude. Tak- 
ing into account the interaction between the mentioned 
tunneling systems gives the right order of the magnetic 
field. 

It should be cited that due to the new experi- 
ments on glasses with different isotopes the effect of 
the nuclear quadrupole moment (NQM) is of significant 
importance 12 which is lacking in our model. The the- 
ory which is based on NQM works with a single molecule 
as a tunneling center and takes its NQM into account. 
However, the next step is to introduce an effective nu- 
clear quadrupole moment for a tunneling system which 
is a cluster of molecules. The theory explains very well 
the low and hight magnetic field behavior of the echo re- 
sponse. But the intermediate regime is missing. Besides, 
in this approach the static electric susceptibility is six 
orders of magnitude smaller than the experimental one. 
We should mention that our model gives the tempera- 
ture dependence of the echo amplitude in the presence 
of a magnetic field which is absent in the NQM model 
proposed in Ref. 0] . 
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By these assumptions, the amplitude of the two pulse 
echo from an effective TLS at time 2t§ measured from 
the first pulse is given by 

A(2t Q ) = (<p(2t ) | p | ^(2i )) 

~ 2Re[a* 1 {T 2 )a 2 {T 2 )p 12 e- 1 ^ SL }, (A2) 

where p = poa z and t 2 = 2t\. 
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APPENDIX A: THE ECHO AMPLITUDE AFTER 
THE SECOND PULSE 

As mentioned in the text the Boltzmann weight is used 
to define the initial condition, |eti(0)| 2 = j?e~ l3El and 
|a 2 (0)| 2 = |e^ B2 , where Z = e"^ 1 + e"^ 2 is the 
partition function, f3 = K T and Kb is the Boltzmann 
constant. So the wave function after the relaxation time 
t is 

| 4>(to)) = ai(n)e h | tpi) + 

-iE 2 (t +T 1 ) 

a 2 (Ti)e s i^). 

During the second pulse the probability amplitudes 
must satisfy Eq. ©, so in a time t' after the second 
pulse we have 

I V^o + t')) = a 1 (r 2 )e^ 1 ^ | <^i) + 
a 2 {r 2 )e ~ \ ip 2 ), 



APPENDIX B: THE DISTRIBUTION OF 
TUNNELING PARAMETERS IN COUPLED 
INTERACTING ETLSS 

If we consider that the interaction between two ETLSs 
is the same as the interaction between the original TLSs 
then we can use the same procedure used previously for 
a pair of TLSs. From Eq.JSJ we have 

fdA 02 f dE 2 E 2 

j J'BS(A r - \Ei- E 2 |H(A„ P - ||^|f ), 

where n(Ei) = [1 +exp((3Ei))~ 1 is the probability of find- 
ing an ETLS in the exited state. By a little manipulation 
we can show that 

T 9T 

f( 3 >(A p ,A 0p )(x-^[ln(- )] 2 . (Bl) 

Op ^Omin 

Since Ao m m ~ T, we find that the distribution function 
for the coupled ETLS parameters is similar to the second 
distribution function (Eq.©). 
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